Surface imaging of inelastic Friedel oscillations 
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Impurities that are present on the surface of a metal often have internal degrees of freedom. 
Inelastic scattering due to impurities can be revealed by observing local features seen in the tunneling 
current with scanning tunneling microscope (STM). We consider localized vibrational modes coupled 
to the electronic structure of a surface. We argue that vibrational modes of impurities produce Fermi 
momentum kF oscillations in second derivative of current with respect to voltage d 2 1(r,V)/dV 2 . 
These oscillations are similar to the well known Friedel oscillations of screening charge on the surface. 
We propose to measure inelastic scattering generated by the presence of the vibrational modes with 
STM by imaging the d 2 I/dV 2 oscillations on the metal surface. 

PACS numbers: 68.37.-d, 72.10.Di, 68.37.Ef 



I. INTRODUCTION 

Inelastic electron tunneling spectroscopy (IETS) with 
scanning tunneling microscope, IETS-STM, is a well- 
established technique, starting with an important exper- 
iment by Stipe et alX These experiments show step like 
features in the tunneling current and local density of 
states (DOS). The physical explanation of the effect is 
straightforward: once the energy of the tunneling elec- 
tron exceeds the energy required to excite local vibra- 
tions, there is a new scattering process that contributes 
to the scattering of electrons due to inelastic excitation of 
the local mode^ The local vibrational modes can be seen 
in differential conductance measurements as side-bands 
to the main (elastic) conductance peak* 3 - Recently, IETS- 
STM was used to measure the spin excitations of indi- 
vidual magnetic atoms4 In addition, recent experiments 
show dominant inelastic channels which are strongly spa- 
tially localized to particular regions of a molecule. 5 

Improvements in spectroscopic and microscopic mea- 
surements have provided new information about funda- 
mental aspects of scattering and interactions in the solid 
state. Recently, STM technique has evolved from imag- 
ing of surfaces^ to spin polarized tips for magnetic sen- 
sitivity of the read-out 4 and spin-polarized injection of 
the tunneling current used in the measurements £ STM 
has been also used to detect Kondo interactions between 
conduction electrons and single atomic spins^ and to 
study the properties of individual atoms i 10 i 1:L 

We propose to use resolution of STM to address spa- 
tially resolved inelastic tunneling features produced by 
local inelastic scattering of the molecule both near the 
scattering center and at far distances. The "fingerprint" 
of the inelastic scattering will be present even away from 
the scattering center. Using local IETS-STM would en- 
able measurements of the standing waves produced by in- 
elastic scattering off impurities, which could be revealed 
as waves seen in the second derivative d 2 I{r,V)/dV 2 . 
For the simplest model of parabolic conduction band 
these waves will be seen as standing waves with period 



set by the Fermi momentum hp. These standing waves 
are seen in the oscillations of the d 2 I(r, V)/dV 2 and are 
similar but qualitatively different from the conventional 
Friedel oscillations. 

In case of regular Friedel oscillations the charge that 
screens off the impurity exhibit oscillations at large dis- 
tances from the impurity. Recent STM measurements 
have observed Friedel oscillations in the elastic scatter- 
ing channels on surfaces with atomic impurities adsorbed 
on the surface These oscillations are seen in a wide 
range of bias as they reflect screening of the charge by 
electron states in the whole bandwidth. Inelastic scatter- 
ing of surface electrons off the molecules may be viewed 
as inelastic Friedel oscillations produced by the electron 
states that are involved in screening. Inelastic Friedel os- 
cillations are seen only in the narrow window of energies 
near the energy of the mode ivo at which inelastic scatter- 
ing occurs, in contrast to conventional Friedel screening. 

It was recently proposed to use IETS-STM to ad- 
dress the inelastic scattering features in the Bi2212 
superconductors i 13 ' 14 In the case of superconductors, the 
situation is more involved since one has to deal with 
the more complicated band structure of Bi2212 and with 
the inelastic scattering off the distributed bosonic modes. 
The work presented here allows one to test similar ideas 
in a more controlled set up where the metal surface and 
molecular modes are well understood. 

The paper is organized as follows. In Sec. fll] we in- 
troduce a model which describe the addressed issue and 
some numerical results are discussed in Sec. IIII1 The 
paper is concluded in Sec. [IV] 



II. PROBING THE FRIEDEL OSCILLATIONS 
AND INELASTIC SCATTERING 
MEASUREMENTS 

The system we consider consists of a two-dimensional 
surface on which inelastic scattering centers are randomly 
distributed. For simplicity, we assume that the impuri- 
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ties are distributed sufficiently far from each other so 
that their mutual interactions can be neglected. All vi- 
brational modes have energy loq and are assumed to be 
the same as they would come from the same type of 
molecules. We use the Hamiltonian for the local vibra- 
tional modes, coupled to electrons via Holstein coupling 
with interactions assumed to occur only at the single im- 
purity site, so that 

U = ^ £kcLr C k<r + (JQ^b + \Y1 C k<T C k',r(^ + b). (1) 



k<7 



kk'cr 



Here, a surface electron is created (annihilated) by 
c ic<T ( c kcr) at the energy e^. The strength of the electron- 
phonon interaction is given by the parameter A, whereas 
luq is the mode of the bare phonon which is created (an- 
nihilated) by b' (b). 

The features we are considering should be seen in the 
second derivative of the tunneling current with respect 
to the bias voltage V in real space, i.e. d 2 I(r, V)/dV 2 . 
This quantity is directly proportional to the frequency 
derivative of the local DOS. In second order perturbation 
theory, this amounts to taking the frequency derivative 
of the correction to the density of states, SN(r, u), due to 
the influence of the impurity scattering. The real space 
electron Green function (GF) is given by 

G(r, r'; u) = G (r, r'; u) + G (r, 0; w)E(w)Go(0, r'; «), 

(2) 

with the zero GF in two spatial dimensions 



Go(r,r'; W ) = ^G (k;a,K k -( r - r ') 

k 

uj + D 



2n ¥ J (k F \r-r'\)(\og 



uj-D 



(3) 



Here, the bare electron Fourier transformed GF 
Go(k, to) = \/{iuj — £k). Using standard procedure for 
frequency summation, the self-energy E(w) is found as 



E(iw) 



k , rn 



^G (k,iuj + in m )Do(iClm) (4) 



^ 2 E 



n B (uo) + /(£k) . n B (w ) + 1 - /(£k) 



ILU + Ll> - £k 



ILU — LJQ — £k 



where Dq(iu>) = 2loq/(ui 2 — Uq) is the bare phonon GFs. 

Due to these definitions, it is clear that the correction 
to the local DOS is given by 

SN(r,Lu) = -IniG o (r,0,w)I](w)Go(0,r,w) (5) 

Since the effects from inelastic scattering are included in 
the self-energy, we evaluate d8N{uj)/duj, which is directly 
proportional to d 2 I(r, V)/dV 2 . 

Sharp feature in self-energy ImE(w) is connected 
to the inelastic scattering feature. Noting that 
ReGo(r, r', u>) w 0, our calculations simplify to 



5N{t,uj) w (l/7r)G (r,0;cj)[ImE(w)]G (0,r;w). We 
thus find that 

SN(r, to) « -t: 3 J 2 (k F r)ImZ(Lu). (6) 

Analytical continuation, e.g. iui — > u> + i0 + , of the self- 
energy in Eq. ((4]) gives the imaginary part 



ImS(w) = -2tHA 2 — [2n B (LO ) + l 

+f(uj + uj Q ) -f(w-w )]. 



(7) 



The correction to the local DOS provides the oscillation 
in real space. The electron-phonon interaction gives rise 
to an increased local DOS for frequencies \u>\ > loq. 
We find that <9 2 /(r, V)/dV 2 is proportional to 

j^SNfaw) « -^ 3 (|?) Jo(k F r) ^ImE^), (8) 



where 



^lmE(w) = 2tt 2 A 2 ^^ fcosh" 2 + cj ) 



d_ 

doj 



h 4 



■ COSh" 2 UJ )). (!)) 



taking (/3/4) cos" 2 0(w ± w )/2 -> 6(uj ± w ), as T -> 0, 
we expect to find sharp features in d 2 I(r, V)/dV 2 around 
a; = ±wo for low temperatures. Hence, our simplified an- 
alytical calculations show that inelastic impurities that 
couple to the surface electrons give rise to Friedel oscil- 
lations in the real space image of d 2 I(r, V)/dV 2 . 

One important point that we also need to address is 
the role of the dephasing in the inelastic scattering pro- 
cess. We are dealing with the inelastic scattering process 
in which the incoming and outgoing electron waves have 
different energy. Energy transferred to/from the inelastic 
center implies the phase change and hence the interfer- 
ence of the incoming and outgoing waves will be affected 
by dephasing caused by scattering. This dephasing pro- 
cess would occur even at T = and hence has nothing 
to do with the thermal scattering. In the limit of small 
boson energy fkoo <C Ep dephasing is small and is not 
going to destroy the interference of the incoming and out- 
going waves. Qualitatively we can estimate the change in 
the phase as a result of inelastic scattering as S(f> = SeSt. 
Here 5e = hu>o is the energy transferred to/from the elec- 
tron and 8t ~ cl/vf is the typical time for the collision 
of electron with the impurity of size a, assumed to be 
on the unit cell size. Then we estimate a — 2n/kF and 
obtain 



2-kTvjJq/Ef < 2tt. 



(10) 



Therefore the interference between incoming and outgo- 
ing waves will survive as long as energy transferred in the 
collision is small. 
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Figure 1: (Color online) Spatial dependence of <5iV(r — ro, u> = 
wo), for T = 1 K. 



III. RESULTS 

Although the analytical calculations in Sec. QTJ clearly 
shows oscillatory response of the inelastic scattering, we 
provide in this section some numerical results to empha- 
size our proposal. The numerical calculations of 5N(r, u) 
are based on the full expression presented in Appendix 

m 

From Eq. ([8]) it is clear that the oscillations only have a 
radial component, because the assumed electron-phonon 
coupling is rotationally invariant. This feature is mani- 
fested in Fig. [TJ where the upper panel shows the spa- 
tial dependence of SN(r — r ,^o) at the resonance en- 
ergy lo = ujq. The lower panel shows the intensity of 
6N(x — a:o,0; a->o), on the line r — ro = x — xq, which 
corresponds to the dashed line in the upper panel. The 
features in SN(r — r ,o;) depend on the spatial position 
and on the frequency. In Fig. [5J which shows a con- 
tour plot of SN(x — xq,0; ui), we illustrate the fact that 
the inelastic Friedel oscillations occur in a narrow interval 
around to — ±luq. The plot clearly illustrates that there is 
hardly any intensity at all, for frequencies off the phonon 
mode loq. The plot in Fig. [3] illustrates the IETS-STM 
resonance signal on a surface with four inelastic scatter- 
ing centers, analogous to the plot in Fig. [TJ As expected 
from the wave character of the signal, the pattern on the 
surface is because of interfering Friedel oscillations from 
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Figure 2: Spatial and frequency dependence of 
\d8N{r,u)/du\ at T= 5 K. 



the four scattering centers. 

Next suppose that the position of the STM tip is kept 
fixed in space, i.e. letting r — ro = where ro is the 
position of the inelastic scattering center. Then, as show 
in the previous section, the derivative d8N(0,uj)/du) is 
expected to peak or dip at u> = ±ojo, c.f. Eq. ([5]) and 
©, as the frequency is varied. This is illustrated in Fig. 
[4] (solid). These local extrema correspond to steps in 
6N(0,uj), which is also clear in Fig. [U (dashed). 

From Eq. ^ we find that intensity of the peaks in 
\d8N/duj\ is highly temperature dependent, which is il- 
lustrated in Fig. O for two temperatures. The plots 
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Figure 3: (Color online) Calculated IEST-STM resonance, 
e.g. at phonon mode resonance, signal on surface with four 
inelastic impurities embedded. 
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Figure 4: (Color online) dSN(0, lo)/8uj (solid) and 5N(0,uj) 
(dashed) for ujq — 3 meV and T = 1 K (m = h = 1). 
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Figure 5: (Color online) The derivative (\d5N/du>\) for dif- 
ferent temperatures T £ {0.01,0.1} K {solid, dashed}. The 
plots have been artificially broadened. 

clearly shows that the peaks are becomes sharper for 
decreasing temperatures, illustrating the S function like 
behavior at low temperatures. 



IV. CONCLUSIONS 

We propose a scanning technique that allows one to 
visualize the oscillations in the inelastic scattering pro- 
duced by the local modes. Standard way to image IETS 
of the molecule is to measure d 2 I(r, V)/dV 2 locally near 
the scattering center We point out that the "finger- 
print" of the inelastic scattering will be present even away 
from the scattering center. The incoming electron wave 
donates/ absorbs energy from the inelastic scattering cen- 
ter. In the process of inelastic scattering the incoming 
and outgoing electron waves interfere. This interference 
will can be clearly seen in the real space oscillations of 
the second derivative d 2 I(r, V)/dV 2 with characteristic 
momentum /cp. The important point is that for the flat 
density of states the interference between states that un- 
dergone inelastic scattering will only be present at the 
bias eV = TvjJq and not at other energies. This makes 
the proposed effect very different from the conventional 
Friedel oscillations that occur for elastic scattering and 
are present for the wide range of energies as the whole 
band of electrons participates in screening. We also ar- 
gue that as long as the energy difference between these 
states is small compared with the Fermi energy of elec- 
trons, the dephasing in the process of scattering is small 
and hence interference is preserved. 

The technique, proposed here, can be applied to the 
inelastic scattering of the electrons from the local vibra- 
tional modes of the molecules on the surface^ and would 
offer exciting possibilities to image the interference and 
interactions between inelastic scattering processes in the 
ensemble of few scattering centers. In practice the tun- 
neling conductance in real materials will likely have peaks 
in d 2 I(r, V)/dV 2 as a function of position due to band 
structure. One would need to be able to separate these 
peaks from the peaks that arise due to inelastic scatter- 
ing. 
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Appendix A: DERIVATION OF 5N 

In the above analysis we neglected the real part of the 
unperturbed Green's functions. However, as we show 
here, taking this part into account does not significantly 
change the picture. We have 



5N(r,u) 



1 /27rmV 



Jl{k F r){ 7r 2 ImE(w) - 



ImE(w) log 



uj + D 



uj-D 

Since we are interested in the frequency derivative of this expression we find that 



- 27rReE(w) 



log 



D 



D 



d_ 

duo 



5N(v,uj) = -^^^(^{Tr^ImSM-^ImSC^log 



uj + D 



AD 

+Im£M — ^log 



2_ D 2 



uj + D 



uj-D 



2n(j^BeZ(u>) ) Los 



uj-D 
uj + D 



uj-D 



2D 1 
- ' -r—fyl) 



UJ* 



Here, 



ReS(w) = 2ttA 5 



n 2 



= 2tt\ 2 — 1 n B (uj ) log 



n B (uj ) + f(e) n B (uj Q ) + 1 - /(e) 
uj + ujq — e uj — ujq — e 

uj + uj + D 



de 



uj + ujq — D 



n B (uj ) + 1] log 



uj — ujq + D 



uj — ujq — D 



(uj - e) 2 - uj£ 



de 



which gives the derivative 
d 



— ReS(w) = -2^- 



[uj + ujQ-e) 2 (uj-uiQ-e) 2 
n B (uj) t n n B (uj a ) + l 



de 



+ 2uj q J f(e) 



D- 



uj — e 



h 2 \ (uj + oj q ) 2 - D 2 (uj - ujq) 2 - D 2 

I 



[(uj-eY-uj 2 ] 2 



de 



However, since the logarithm log | (uj ± uj + D)/ (uj ± uj — first are negligible in d5N/duj . 
D)\ k for D > \uj + uj \ we see that all terms but the 
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